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ABSTRACT 

We present line strengths in the bulges and inner disks of 38 galaxies in the local 
universe, including several galaxies whose bulges were previously identified as being 
disk-like in their colors or kinematics, to see if their spectral properties reveal evi¬ 
dence for secular evolution. We find that red bulges of all Hubble types are similar to 
luminous ellipticals in their central stellar populations. They have large luminosity- 
weighted ages, metallicities, and a/Fe ratios. Blue bulges can be separated into a 
metal-poor class that is restricted to late-types with small velocity dispersion and 
a young, metal-rich class that includes all Hubble types and velocity dispersions. 
Luminosity-weighted metallicities and ajFe ratios are sensitive to central velocity 
dispersion and maximum disk rotational velocity. Red bulges and ellipticals follow the 
same scaling relations. We see differences in some scaling relations between blue and 
red bulges and between bulges of barred and unbarred galaxies. Most bulges have 
decreasing metallicity with increasing radius; galaxies with larger central metallicities 
have steeper gradients. Where positive age gradients (with the central regions being 
younger) are present, they are invariably in barred galaxies. The metallicities of bulges 
are correlated with those of their disks. While this and the differences between barred 
and unbarred galaxies suggest that secular evolution cannot be ignored, our results 
are generally consistent with the hypothesis that mergers have been the dominant 
mechanism responsible for bulge formation. 
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1 INTRODUCTION 

Bulges are important relics of the galaxy formation process. 
An analysis of their structure, kinematics, dynamics, and 
stellar content can potentially reveal the physical mecha¬ 
nisms responsible for the formation and evolution of galax¬ 
ies as well as the nature of the Hubble sequence. Similarities 
between bulges and ellipticals have long been recognized but 
recent observations suggest that at least some bulges may 
be related to disks. This has led to the suggestion that the 
large bulges of early-type spirals are more sim ilar to ellip¬ 
ticals while late-type bulges are more disk-like JWvse et alJ 
Il99it) . As a consequence of these observations, formation 
scenarios have emerged for bulges that are either identical 
to those for ellipticals or involve the secular evolution of 
disks. However, the degree to which formation mechanisms 
are homogeneous is still open to question. 


* Based on observations obtained with the Apache Point Obser¬ 
vatory 3.5-meter telescope, which is owned and operated by the 
Astrophysical Research Consortium, 
f E-mail: bmoorthy@nmsu.edu 


Early models for elliptical formation invol ved the mono¬ 
lithic collapse of a primordial gas c loud llLarsorJ Il974l : 
ICarlber3ll98dL lArimoto fc Yoshiiill987ll . This model natu¬ 
rally explains several observed properties of ellipticals in¬ 
cluding the mass-metallicity relation and the presence of 
metallicity gradients but large-scale collapse is inconsis¬ 
tent with present day cold dark matter cosmology and 
with recent observations showing that m assive ellipticals 
were not fully as sembled since after z=l llBell et al.ll2004l : 
iFaber eTaOEo^. It is now widely believed that ellipticals 
formed hierarchically th rough mergers of smaller fragments 
iKaiiffrn ann et ^Il993|). Mergers are frequently caught in 

the act Tvm^DokMrii^t alll99fl( : lFerreiro fc PastorizalAo^l 

and photometric and ki nematic evidence for past mergers is 
abun dant in ellipticals iEmsellem et aljl2004l : Ivan DokkumI 
[ 200 ^. The merger model has been extended to bulges due 
to the many observ ed similarities betwe en bulges and ellip¬ 
ticals. For example lCarollo et al.l lll997^ found that bulges 
were well-fit by the law used for ellipticals. The funda¬ 
mental plane relation of bulges is nearly the same as that of 
ellipticals, with late-types p erhaps lying below early types 
iFalcon-Barroso et alJiAo2ri . 
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In the secular evolution scenario, bulges are produced 
through radial and vertical transport of disk material as 
the result of instabilities and resonances (see Kormendy & 
Kennicutt 2004 for a review). These models come in sev¬ 
eral flavors, most of which involve bars. Simulations that 
do not include gas have found that bars can buckle, heat¬ 
ing the inner disk and increasing its scale height to resem¬ 
ble a bulge. Hydrodynamical simulations have found that 
bars can transport g as towards the center, triggering in- 
tense star-formation JPfeimigeijJlOOjjFri^^ 

iNorman et al ]ll990l:lNoguchil200nl:llmmeli et alJl2004il. The 

S nce of neighbors may also drive this iKannappanet alJ 
. Support for secular evolution comes from observed 
correla tions between the sca l e lengths of bulges and their 
disks iCourteau et alJ 119961 : iMacArthur et all l200,?ll . Re¬ 
cent work has also shown that the light pr ohles of many 
bulges are closer to exponential than Balcells et alJ 


3 3S are closer to exponent i al tuan K ' cmce^ 
; IMacArthur et all 1200.^ : Ide Jong et al ~ l20(^r 


Fur¬ 
thermore, the ratio of rotationa l to random motions in 

S 5S is often typical of disk s iKor mendv fc IllingworthI 
; iKormendv fc Kennicutli l20o3h Comparisons be- 
tween the morphology and kinematics of observed galax¬ 
ies with simulated ones have shown that boxy and 


clination JBureau & FreemarJ Il999l: 

Aronica et alJ 

l200,3l: 

iGhune & Buread l2004l: lAthanassoulal 

200,. Athanassoula 


(2005) distinguishes between b/p bulges, which are formed 
through the buckling of the bar, and what she calls “disky 
bulges”, which are smaller cold components that formed out 
of the gas driven inward by the bar. 

Bulges that could have been formed through secular 
evolution are often referred to as “pseudobulges” to distin¬ 
guish them from the “classical” bulges that may have formed 
through mergers. Since pseudobulge signatures are generally 
found in later-typed spirals, Kormendy & Kennicutt (2004) 
suggest that early-type spirals (Sa’s, Sab’s, and some SO’s) 
contain classical bulges while late-type spirals (Sb’s, Sc’s and 
some SO’s ) cont ain pseudobulges. On the theoretical side, 
IPfennig^ lll99.?ll found that secular evolution can produce 
small bulges but not those having a characteristic radius 
much larger than the disk scale length. However, it is not 
at all clear that the spectrum of observed bulge properties 
points towards two distinct formation sc enarios. Since the 
stability of bars continues to be debated _l |Shen_fc SellwoodI 

l2004l : lDebattista et al.l2004l : lBournaud eTahl^OsIl . it is also 

not clear whether or not pseudobulges should exist only in 
present-day barred galaxies. 

Stellar population (SP) studies can potentially place im¬ 
portant constraints on the formation mechanisms. A succes- 
ful formation scenario has to reproduce the observed distri¬ 
bution of ages and metallicities. In a collapse model, bulges 
and ellipticals are universally old and have radial metallic- 
ity gra dients. In his dissipative collapse simulation ICarlberd 
^1984^ found that the steepness of the metallicity gradient 
was correlated with galaxy properties such as mass and lu¬ 
minosity. If ellipticals and bulges formed through mergers, 
it is important to keep in mind that their assembly histo¬ 
ries might be very different from their star formation his¬ 
tories. ACDM simulations suggest that most massive ellip¬ 
ticals (and therefore presumably bulges) were not fully as¬ 
sembled until recently (z<l) whereas the bulk of star forma¬ 
tion occurred much earlier (z>2) in the progenitor galaxies 


jPe Lucia et al 11^ - This is consistent with observational 
studies of merg er activity, numbe r counts, and the l umi- 
nosity function llFaber et al.ll200^ Ivan DokkumlEoO^ . de 
Lucia et al. find that the star formation histories of massive 
ellipticals peak at z«5 while those of less massive ellipti¬ 
cals peak at progressively smaller redshifts and are more 
extended. These simulations predict a mass-metallicity rela¬ 
tion, with the most massive ellipticals having solar metallic¬ 
ity and the least massive ones being a factor or ten smaller 
in metallicity. Gradients in SPs are difficult to model within 
the framework of hierarchical formation. Mergers between 
disks can presumably preserve existing gradients in the sub¬ 
components and produce new gradients through gas infall, 
but mixing from successive mergers might erase any correla¬ 
tions between gradients and global properties. White (1980) 
found that the metallicity gradient in a disk galaxy was 
halv ed after three mergers with similar sized disks. How¬ 
ever, ^ekk^^^h^^ lll999^ found that more massive galax¬ 
ies had steeper metallicity gradients. The impact of secu¬ 
lar evolution on gradients is not straightforward. Since the 
resulting pseudobulge has a smaller scale length than the 
progenitor disk, an existing disk gradient could become am¬ 
plified (A Klypin, private communication). However, mixing 
during secular evolution can have the effect of washing out 
existing gradients. Adding gas only complicates the picture. 
If gas is fueled towards the central regions by bars, this could 
result in a nucleus that is younger and more metal-rich than 
the outer regions of the bulge. Simulations by Friedli et al. 
(1994) resulted in a flattening of metallicity gradients in all 
but the innermost regions where a starburst, fueled by in¬ 
falling gas, produced a metal-rich nucleus. 

Abundance ratios can place additional constraints. Mg 
and other a-elements are primarily produced in Type H Su¬ 
pernovae (SN H) while a substantial fraction of the Fe-peak 
elements Fe and Cr are produced in Type la Supernovae 
(SN la). Therefore, a-enhancement is generally attributed 
to a cessation of star formation before the bulk of SN la 
occurred. Through chemical evolution modeling, Thomas et 
al. (1999) found that a clumpy collapse model produced uni¬ 
form a-enhancement or positive gradients (increasing a/Fe 
with radius) while a merger model produced uniformly so¬ 
lar g/ Fe or negative gradie nts. For the case of secular evo¬ 
lution, llmmeli et al.l l^2004^ make different predictions for 
abundance ratios in bulge stars depending on whether it is 
the gas disk or the stellar disk which first becomes unstable. 
In the former case, gas clumps merge together and spiral in¬ 
ward, causing massive starbursts and producing large a/Fe 
ratios. In the latter case, a bar forms and then channels gas 
towards the center. This occurs on long timescales, resulting 
in smaller a/Fe ratios. 

The only bulges where individual stars can be resolved 
are those of the Milky Way (MW) and M31. The major¬ 
ity of the stars in the MW bulge are old (t^7 Gyr), al¬ 
though young (ts;200 Myr) and int ermediate-age (200 Myr 
4t^7 Gyr) stars are also detected ilbat.LJ^TjjImQrj|l9^: 

SadIerjjtjaJjy^99^ J^ fc Gilmorel200nnvar^oone^l] 

20031: IZoccali et alJl2nn3lL As a barred late-type spiral, the 

MW might be a good candidate for secular evolution but 
that hypothesis is challenged by the mean stellar age of 
its bulge. If the bulge were produced through a rearrange¬ 
ment of disk stars, this must have occurred several Gyr 
ago if the inner disk has the same age distribution as the 
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disk at the solar neighborhood. The young stars in the 
bulge are mainly found in the innermost regions. While 
the gas that formed them could have been driven by a 
bar, it could just as well have been provided by a re¬ 
cent merger. The mean metallicity of the Galactic bulge is 
slightly sub-sola r <lMinnitilll996l: IPeltzing fc Gilmor j I 2 OOOI : 

Ivan Loon et al ] l2003l : IZoccaiietahll20n.'jh Recent Fullbright 

et al. 2006). The stellar content of M3Ts bulge is not as well 
understood as that of the MW since we cannot reach its main 
sequence turnoff. Observations of giant stars are consistent 
with M3Ts bulge being si milar in age to the MW bulge and 
slightly more metal-rich iRich ||l99d |Jablonka^.et_alJ I 2 OQOI; 

S iral200lt ISteohens et all 1200?ISa^edhi^^abloni^ 

SPs in more distant bulges have to be studied through 
photometry or spectroscopy of integrated light. An im¬ 
portant limitation of such studies is that integrated light 
is dominated by the most luminous stars. Colors have 
been studied more extensively as they have the advan- 
tage of higher signal-to- noise (S/N). Pioneering work by 
iBalcells fc PeletiCT if994) found that color variations from 
galaxy to galaxy are much larger than color differe nces be¬ 
tween disk and bulge in each galaxy. Similarly, Ide .Tond 
il99fil found that bulge and disk colors are correlated and 
that the color differences between bulge and disk suggested 
that the SPs did not vary much from one to the other. Un¬ 
fortunately, color studies suffer from degeneracies between 
ages, metallicities, and extinction. 

Line strengt hs are nearly insensitive to dust 
iMacArthurll200.^ . provide information on the abundances 
of several elements and molecul es, and a ll ow fo r breaking 
the age-metallicity degeneracy. IWorthe-^ (Il994l obtained 
line strengths for a range of single a ge, single metallicity 
SPs (SSPs) on th e Lick/IDS system IIBurstein et al.lll984 
iFaber et al.lll98!ill and found that while individual indices 
are sensitive to both age and metallicity, the relative 
sensitivity varies from index to index. Spectral indice s 
have also been dehned at high resolution llVazdeki£lll994l . 
allowing better age determinations than otherwise possible. 
One of the limitations of the original models is that they 
were calibrated using galactic stars, few of which had 
abundance ratios different from solar. Much progress has 
since been made in extending Lick indices to non-solar 
abundance ratios_jTVigic^^& Belll99,4:lTrager et al.l2000al: 

iThomas et al 1120031.12004 iLee fc Worthevll200,4l. 

Line strengths have been used extensively to charac¬ 
terize the SPs of ellipticals. The luminosity-weighted ages 
of cluster ellipticals are large while those of field ellipticals 
are on average smaller, with a large sprea d (|Rose_et_aI .11994 

[Tragergtjalj2000.jjVazdekis_gt_alj200i y^roctCT^^SansM 

l200?lDenicolWe^jTl200^lThoma^^lTl200RlL ~ This goes 

against the collapse model and conhrms, at least qual¬ 
itatively, the prediction of the merger model by Kauff- 
mann (1996). Mg-sensitive indices in ellipticals are more 
tightly correlated with central velocity dispersion than Fe- 
sensitive indices, resulting in a correlation between Mg/Fe 
and CTO llBender et alJ Il993ll . IWorthev fc Colloberd (l2003ll 
found that that the Mg-cr relation of ellipticals is consistent 
with these objects having been formed through around 50 
mergers with merger probability constant or midly declining 
with time. 

There have been fewer studies of line strengths in 


bulges. Integrated light studies on the bulges of the MW 
and M31 have ar rived at similar ages and metallicities as the 
resolved studies llPuzia et alJl200^l2^^5^ . Both bulges have 
large SSP ages. M31 is slightly super-solar in SSP metal¬ 
licity while the Milky Way is solar. Both are a-enhanced 
with M31 being more so in line with its larger ctq. Early 
studies on extragalactic bulges found the m to be similar 
to ell i pticals in their cen t ral line streng ths jjablonka e t all 
Il996l : lldiart et al.lll99(fl . IProctor Hr. SansomI (l2002|l found 
that bulges have smaller average luminosity-weighted age 
than ellipticals. These authors did not Hnd the correla¬ 
tion predicted by Kauffman (1996) for bulges and suggested 
that it might have been erased by secular evolution in late- 
types. The largest s ample of bulges to date was that of 
IPrugniel et al.l 112004) . who identihed three classes of bulges: 
a) young bulges which are small, have ionized gas, low ve¬ 
locity dispersions, and low metallicity; b) old bulges that 
are alpha-enhanced and follow the mass-metallicity relation 
of ellipticals; and c) bulges that have a mixture of young 
and old populations, which are less alpha-enhanced than 
those of class (b), and deviate from the Mg 2 relation of el¬ 
lipticals. Prugniel et al. and Proctor et al. found that both 
Fe and Mg were correlated with ctq in bulges, resulting in 
the lac k of a tight co r relatio n between Mg/Fe and ctq in 
bulges. IPrugniel et al.l (1200 ill found that Mg 2 in bulges is 
more tightly correlated with the Ymax of the disk than with 
CT, indicating that the SPs are more sensitive to the total 
galaxy potential (i.e. the dark matter halo) than the bulge 
potential. 

Studies with spatial resolution offer several advantages 
to studies that only sample the central region. First, dif¬ 
ferential studies of ages and abundances are more reliable 
than absolute estimates. Second, formation models invari¬ 
ably make predictions for the global properties of galaxies 
which are better traced by mean observed quantities than 
central ones; observations with spatial resolution allow esti¬ 
mation of mean values. Finally, as mentioned already, pop¬ 
ulation gradients can place additional constraints on forma¬ 
tion mechanisms. 

Line s trength gradients ha ve b een studied ex t ensive ly in 
ellipticals. ICarollo et al.l (ll993^ and iForbes et al.1 (1200,4) hnd 
strong correlations between gradients and physical prop- 
erties while others hnd wea k (Mehlert et al. 2003) or no 
iKobavashi fc Arimotolll999^ correlations. There have been 
relatively few studies on gradients in bulges. Fisher et al. 
(1996) found steeper metallicity gradients along the mi¬ 
nor axes of nine edge-on SOs than along the major axes, 
suggesting di herent formation mech anisms for the bulge 
and the disk. iGoudfrooii et all (Il999l) found that gradients 
were c orrelated with luminosity in 16 bulges. IProctor et alJ 
ll200ri) found that gradients correlated with velocity dis- 
persion, albeit with a sample of only four galaxies, while 
I.Tablonka et alJ (l2002l) found no such correlation. Integral 
held spectroscopy has enabled the acquisition of 2D line 
strengths in bulges with results just starting to emerge (e.g. 
Sil’chenko et al. 200 3; Falc on-Barroso et al. 2004). Recent 
work bv iR vder et al.l (1200,4) shows that tunable hlters might 
be another way to obtain 2D line strengths. 

In this paper, we present line strengths and line strength 
gradients in the bulges and inner disks of 38 galaxies. Our 
sample, described in Section 2, was chosen to span a range 
of bulge properties and specihcally targeted several galaxies 
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with blue bulges and similar bulge/disk colors and/or disk¬ 
like kinematics in an attempt to look for SP signatures of 
secular evolution. Section 3 describes the observations and 
data analysis. Section 4 describes the SP results and Section 
5 discusses their implications for bulge formation scenarios. 
Section 6 contains a summary. The structure, kinematics, 
and dynamics and how they relate to the SPs will be dis¬ 
cussed in a future paper (hereafter Paper II). 


2 THE GALAXY SAMPLE 

We selected a sample that included some bulges that are 
similar in color to their disks and others that are consider¬ 
ably redder as a control. Color was chosen as the primary 
selection criterion becau se this ha s so far been the best stud- 
ied property of bu lges. Ide Jonel lll99fil . hereafter DJ) and 
IPeletier fc Balcella il997l hereafter PB) obtained col or gra¬ 
dient s of galaxies from the Uppsala General Catalog (iNilsonI 
Il97,'jl with major axes larger than 2’. We selected 17 galaxies 
PB and 14 from DJ. The two samples complement each other 
nicely in their sky coverage and sampling of Hubble types. 
The DJ galaxies are nearly face-on while the PB galaxies 
are highly inclined. 

We also included three galaxies, NGCs 2787, 3384, 
and 3945, which were previously fou nd to possess disk-like 
struc tural and kinema ti c properties (|&ugMelIo_et_alJ[l996t 

S enko et ^ l2003t IPinknev et alJ l20oS ~~ lE^ine^l] 
. All three are barred SO galaxies with inner 
disks or bars that are more luminous than the sur¬ 
rounding bulge. One of the PB galaxies, NGC 7457, 
is al so known to h a ve di sk- like kinematics iKormen d3 
1 19931 : IPinknev et alJ 1200311 . iMichard fc MarchaJ 1^941 


found small bar-like distortions in this galaxy and 
lEmsellem et alJ (l2004^ found nearly cylindrical rotation 


which is seen in boxy bulges ( 

Bertola & Canacciol 

197';i 

iKormendv V. IllingworthI Il982l: 

Falcon-Barroso et al 

2004[) 

and in simulations of edge-on 

bars JCombes et al. 

199qI 


ISellwoo^ll993l : lAthanassoula fc Misirioti^l^OO^T ' 


This project initially began in collaboration w ith some 
members of the ENEAR survey llWegner et al.l200(ll . There¬ 
fore, the first hve galaxies we observed were from their sam¬ 
ple: NGGs 4472, 2775, 3544, 3831, and 5793. NGC 4472, a 
bright elliptical in the center of the Virgo cluster, was in¬ 
cluded for comparison with previous studies. The other four 
galaxies were selected to span a wide range in inclination 
and bulge-to-disk ratio. 

There are several reasons for including both high- and 
low-inclination galaxies: 

1. Minor axis observations of highly inclined galaxies 
offer minimum disk contamination in the outer regions of 
the bulge. 

2. In moderately inclined galaxies, there is actually more 
disk contamination along the minor axis than the major 
axis for the same solid angle. To estimate the degree of disk 
contamination, we obtained spectra along both major and 
minor axes for some of our inclined galaxies. 

3. Major axis observations of inclined galaxies allow for 
the measurement of rotation which can provide additional 
information about the structure of the galaxy. 

4. Low-inclination galaxies have less disk contamination 
in the central regions and allow for clear identification of 


bars, rings, and other morphological features. Including both 
high- and low-inclination galaxies allows for a comparison 
between bars and b/p bulges. 

5. In highly-inclined galaxies, the bulge and disk can 
be distinguished based on their shapes (spheroidal versus 
flat). In face-on galaxies, this is not possible and so bulges 
are generally defined as the excess light on top of the inward 
extrapolation of an exponential disk. SPs and kinematics of¬ 
fer two additional and independent means of distinguishing 
between bulges and disks in face-on galaxies. 

Twelve out of our 20 low-inclination galaxies are barred. 
Some of our highly -inclined galaxies were classified by 
iLiitticke et alJ ll2nnfi) into peanut-shaped, boxy, nearly boxy, 
or elliptical bulges. For our remaining highly-inclined galax¬ 
ies, we determined the shapes using their technique. This 
yielded 10 b/p bulges and 8 elliptical bulges. Therefore 
the fraction of barred galaxies in low-inclination galaxies 
is approximately equal to the fraction of b/p bulges in 
highly-inclined galaxies. While peanut-shaped bulges are 
easily identihed, it is not always easy to distinguish an el¬ 
liptical bulge from one that is slightly boxy. For instance, 
Ltitticke et al. class i fy NG C 5838’s bulge as elliptical but 
iMichard fc MarchaJ (ll994^ describe it as boxy. 

Table 1 contains basic data on our galaxies. The column 
“Morph.” describes the shape of the bulge if the galaxy is 
highly-inclined (Boxy, Peanut, or Elliptical) and whether or 
not it is barred if it is not highly-inclined. Identifications 
marked with an asterisk are those of Lfitticke et al. while 
those without asterisks are our identihcations. 

When comparing SPs in galaxies with different colors, 
it is important to keep in mind that color is correlated with 
the global dynamical properties of a galaxy. Fig. 0 shows 
the bulge B-K colors as a function of central velocity disper¬ 
sions, and maximum disk rotational velocities of our galax¬ 
ies where available. We found that it is useful to subdivide 
bulges according to whether they are redder or bluer than 
B-K=4; these are shown as red and blue points. In this and 
subsequent plots, point shape represents the Hubble type: 
circles are SO’s; hexagons are SOa’s and Sa’s; pentagons are 
Sab’s; squares are Sb’s; and triangles are Sbc’s and Sc’s. 
Filled symbols are barred galaxies. Thin open symbols are 
elliptical bulges if highly-inclined and unbarred galaxies if 
not highly-inclined. Thick open symbols are b/p bulges. 

Bulge colors correlate more tightly with Ymax than with 
(To. Galaxies with Vmax > 200kms~^ host red bulges while 
those with Vmax < 165fems“^ host blue bulges. Both red 
and blue bulges are found in nearly the full range of central 
velocity dispersions spanned by our galaxies although there 
is an overabundance of red bulges in large a galaxies and 
vice versa. 

3 OBSERVATIONS AND DATA ANALYSIS 
3.1 Observations 

Observations were made with the Double Imaging Spectro¬ 
graph (DIS) on the ARC 3.5m telescope at Apache Point 
Observatory between January 2000 and February 2004. The 
spectrograph uses a dichroic to split the light into separate 
blue and red channels. During this period, the instrument 
was upgraded in several phases with the installation of new 
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Figure 1. Bulge color versus central velocity dispersion and max¬ 
imum disk rotational velocity. The velocities are from this work 
while the colors are from PB and DJ. The central velocity dis¬ 
persions were measured within an aperture of approximately 4 
arcsec. See Table 1 for the definition of bulge color. Bulges are 
shown in red or blue according to whether they are redder or bluer 
than B-K=4. Filled symbols are barred galaxies. Thin open sym¬ 
bols are elliptical bulges if highly-inclined and unbarred galaxies 
if not highly-inclined. Thick open symbols are b/p bulges. 


detectors and optics. Table 2 gives the specifications of each 
configuration and Table 3 describes the spectroscopic obser¬ 
vations. DIS I gave us continuous wavelength coverage from 
4000 to 7500 A while DIS II and DIS III gave us continuous 
coverage from 3700 to 7500 A. A 5 arcmin x 1.5 arcsec slit 
was used in all the observations. On each night, we observed 
a quartz lamp for fiatfielding, arc lamps for wavelength cal¬ 
ibration, and two to five spectrophotometric standards for 
flux calibration. On s everal nights, Lick standard stars from 
IWorthev et~ai1 J 19941) were observed to allow us to transform 
our line indices to the Lick/IDS system. 

For most of the highly inclined galaxies, we obtained 
spectra along both major and minor axes. For the unbarred 
low-inclination galaxies, we obtained major axis profiles ex¬ 
cept for NGC 2916 for which we obtained a minor axis proHle 
instead because there was a bright star along the major axis. 
For the clearly barred galaxies, we placed the slit along the 
bar. For IC 302 and NGC 2487, this is different from the 
position angle of the major axis. In NGC 5375, the bar hap¬ 
pens to be along the major axis. We could not identify a 
major axis in NGC 266 or NGC 5020 but the slit must have 
been placed close to it as we see substantial rotation. We did 
not detect any rotation in five galaxies: IC 302 and NGCs 
765, 2487, 2916, and 6246A, due to their inclinations being 
too low. 

We obtained images in B, V, and R for bulge-to-disk 
decomposition on six nights using the SPIcam detector on 
the same telescope. Typical exposure times were 600 sec for 
B, 180 sec for V, and 120 sec for R. Superbiases and twilight 
flats were obtained on each night. 


3.2 Basic Reductions 

Data reduction was carried out using the XVISTA soft¬ 
ware package. For the imaging, basic reduction included bias 
subtraction and fiatfielding. For the spectroscopy, flat-helds 
were constructed using a median of 5 to 10 bright quartz 
lamp exposures; the mean spectral response was divided out. 
Wavelength calibration was performed using Fie, Ne, and 
Ar arc lamp exposures, using a hfth order polynomial for 
both blue and red channels. Flux calibration was performed 
using a spline fit to published spectr a of the spectrophoto¬ 
metric standards llMassev et alJll988^ . Line curvature along 
the slit was measured using the lamp exposures and a sim¬ 
ple row by row shift was stored for subsequent correction. 
Similarly, spatial distortion in the spectrograph was mea¬ 
sured using standard star exposures and the correction was 
stored. Since spatial distortion includes a component due 
to differential refraction unless the slit is perpendicular to 
the parallactic angle, this component was calculated and re¬ 
moved from the standard star measurements. Application of 
a correction to the galaxies includes a refraction component 
and the spectrograph component. 

Each galaxy frame was subtracted by the overscan and 
superbias, flat-helded, and trimmed to remove the overscan 
region. The line curvature and spatial distortion corrections 
were applied. Multiple observations of each galaxy were then 
coadded, rejecting cosmic ray outliers in the process. Where 
multiple observations were not available, cosmic rays were 
removed by eye using a spatial median Hlter. Variance frames 
were propogated throughout the reduction process. To com¬ 
bine the red and blue channels, the red galaxy spectra were 
rescaled to spatially match the blue spectra. The spectra 
were then extracted in approximately 1 arcsec bins near the 
center of the galaxy and using larger bins further out. An av¬ 
erage of sky values measured on both sides of the galaxy was 
subtracted from each slice. Wavelength and flux calibrations 
were applied and both red and blue frames were rebinned 
to 3 A/pixel. Finally, the spectra were deredshifted and the 
red and blue sides were combined. The blue spectrum was 
used out to a wavelength of 5600 A and the red from 5450 
A. In the overlap region between 5450 and 5600 A, an aver¬ 
age of red and blue values was used. In some galaxies, there 
is a discontinuity in the overlap region that arises from the 
difficulty of accurately combining the two channels for ex¬ 
tended objects; fortunately none of our absorption features 
fall within this region. 

On our last observing run (Feb 15, 2004), the light 
on star frames was too extended to be explained by seeing 
alone. The excess light, which we believe is scattered light, 
had a spectral energy distribution similar to that of the 
mashed stellar spectrum but without any narrow spectral 
features presumably because it is significantly defocused. 
Since adding a constant to a spectrum decreases the equiv¬ 
alent width (EW) of an absorption features and since the 
relative contribution of the scattered light to galaxy light in¬ 
creases with distance from the galaxy center, this introduces 
an artificial negative gradient in the line strength proHles. 
To correct for the scattered light, we fit the 2d stellar spec¬ 
trum with a smoothed stellar spectrum along the wavelength 
direction and a fifth order polynomial along the spatial di¬ 
rection, masking out the central 20 pixels (8.4 arcsec). This 
spatial profile, combined with the smoothed spectral prohle 
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of the galaxy, was subtracted from each galaxy frame. One 
galaxy that was observed on the problematic night, NGC 
3384, has previously measured index prohles. Applying the 
correction resulted in much better agreement with the pub¬ 
lished values (see Section 3.7). No scattered light correction 
was applied on any of the other nights since the correction 
derived for those nights did not affect the line strength pro¬ 
files significantly. 

3.3 Measuring and Correcting for Rotation and 
Velocity Dispersion 

We measured rotation and velocity dispersion in the stel¬ 
lar components of our galaxies using the pPXf pack¬ 
age llCa.nnellari V Em.sel1eml l2004^. SSP spectra were con- 
structed using the SP models bv iBruzual fc Charlotl i2003L 
hereafter BC03) assuming a Chabrier initial mass function. 
The pPXF routine fit each galactic extraction with a lin¬ 
ear combination of SSP spectra, shifting and broadening 
these to match the galaxy’s rotation and velocity dispersion 
respectively. The fit was performed within the wavelength 
range 4800-5400 A, with the emission lines H/3, [OIII] 4959, 
and [OIII] 5007 masked out. The profiles will be presented 
in Paper II. 

The galaxy spectrum at each location was shifted by the 
measured stellar rotation before measuring the line indices. 
To correct the indices for velocity dispersion, line strengths 
were measured on an SSP template that was broadened by 
the measured velocity dispersion and an unbroadened but 
otherwise identical template. These templates were also con¬ 
structed using a linear combination of BC03 SSP spectra 
with emission lines masked out, but the wavelength range 
for the fit was 4000-6600 A, to include all the Lick indices. 
This template was also used for emission correction (follow¬ 
ing section). The measured absorption line equivalent widths 
were multiplied by the ratio of the unbroadened line strength 
to the broadened one; for magnitudes the correction factor 
is the difference between these two quantities. 

3.4 Measuring and Correcting for Emission 

Some absorption indices can be severely affected by line¬ 
filling by emission. These include the Balmer indices, Fe5015 
(due to [OIII] 5007 emission), and Mg;, (due to [NI] 5199 
emission). To correct for this, we subtracted the template 
described in the previous section from each galaxy spec¬ 
trum. If on the residual spectrum, Ha was found to be in 
emission at the 5a level and a local maximum was detected 
at H/3, a gaussian was fit to the H/3 emission and subracted 
from the galaxy spectrum. This procedure was repeated for 
Hy and H5. [OIII] 5007 and [NI] 5199 were subtracted out 
if they were found to be in emission at the 3a and ia lev¬ 
els respectively. A larger threshold was used for [NI] 5199 
because this feature lies at the edge of the Mg 2 absorption 
feature and spurious discontinuities often show up there due 
to template mismatch. 

EWs were measured for several emission-lines to study 
the nature of the ionized gas in bulges and inner disks. The 
galaxy continuum, obtained by smoothing the galaxy spec¬ 
trum, was added to the emission spectrum described above 
before measuring the EWs. 


3.5 Lick Index Measurements 

The final galaxy spectra were broadened to 9.5 A FWHM, 
which is approximately equal to the Lick resolution, and 
rebinned to a dispe rsion 0.125 A/pixel. Variabl e broaden¬ 
ing as prescribed in IWorthev fc Ottavianil ^1997^ was tried 
and found not to produce significantly different results. The 
strengths of 25 absorption features were measured using 
the latest bandpasses from Guy Worthey’s webpage. The 
EW or magnitude of each feature was computed following 
iTVa.ger ~ Jl998h . 

Spectral indices were measured on Lick standard stars, 
exactly as done for the galaxies, to transform our line 
strengths to the Lick/IDS system. Such transformations are 
necessary because our detector differs in resolution from the 
IDS and because the IDS spectra were not flux calibrated. 
24 stars ranging in spectral type from F5V to K7HI were 
used inderiving the transformations for DIS I and 22 stars 
(also F5V to K7HI) were used for DIS HI. DIS HI transfor¬ 
mations were used for DIS H since the same gratings were 
used for both. From imaging one of the stars at several posi¬ 
tions along the slit, it was determined that the line strengths 
do not vary significantly with slit position. Fig.|21shows the 
transformations. The rms scatter when applying the trans¬ 
formations on the stars are shown at the bottom-right of 
each plot. The flux response of the detector changes rapidly 
at the long wavelength end of the blue channel, making mea¬ 
surements there difficult. Only the Mgi and Mg 2 indices 
were affected by this, but severely so, since their continuum 
bandpasses are far apart and the red continuum bandpass 
lies right where our flux response is steepest. This resulted in 
large scatter in the transformations of these indices. For the 
other indices, the scatter in the transfo rmations is similar 
or slig htly larger than that obtained bv IProctor fc Samsoml 

J200a. 


3.6 Comparison with the literature 

Fig. 0 shows comparisons of Lick indices between this work 
and the most recently published values for NGG 4472. There 
is excellent agreement in the Ga4227, H/3, Fe5015, (Fe),and 
Fe5270 profiles between this work and all the published val¬ 
ues. Our other index prohles are systematically offset from at 
least one of the other two studies but the offset is within the 
uncertainties of transforming our data to the Lick system. In 
general, our results ag ree better with the long -slit data from 
IVagdekis et alJ lll997^ than the IFU data from IPeletier et ^ 
^1999^ . Nearly all our central values agree with those of 
Trager et al. Emission correction was not responsible for 
any disagreement among studies since NGG 4472 does not 
have much emission. 

The hrst column of Fig. |1] shows index comparisons 
between this work (triangles), azimuthally averaged IFU 
data from Sil’chenko (1999; asterisks), and SAURON data 
extracted along the major-axis from Falcon-Barroso et al. 
(2004; squares) for NGG 7332. All the central values are 
in agreement. The slope of our H/3 and Fe5270 prohles are 
steeper than Falcon-Barroso et al.’s but not as steep as 
Silchenko’s. Our Mg;, and Fe5015 prohles agree reasonably 
well with Falcon-Barroso’s but Silchenko’s Mgb prohle is 
again steeper. 

The second column of Fig. 2] shows index comparisons 
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Figure 2. Transformations to the Lick/IDS system. Each point is a Lick standard star. Triangles are DIS I. Crosses are DIS III. DIS III 
transformation were used for DIS II. The rms scatter when applying the transformations on the stars are shown at the bottom-right of 
each plot (with DIS I above DIS III). 


between this work, long-slit data from Fisher et al. (1996; 
crosses), IFU data from Sil’chenko et al. (2003; asterisks), 
and SAURON data from de Zeeuw et al. (2002; squares) 
for NGC 3384. Here also, Sil’chenko’s profile is averaged az- 
imuthally while de Zeeuw et al.’s is an extraction along the 
major axis. This galaxy was observed on the night in which 
scattered light was corrected for as described in Section 3.2. 
Fisher et al.’s profiles are in good agreement with those from 
SAURON. Our Mgt and Fe5270 profiles agree marginally 
with published values. In the central regions, our H/3 profile 


falls in between de Zeeuw et al.’s and Sil’chenko et al.’s. Out¬ 
side about 15 arcsec, our profile falls steeply, possibly due 
to inadequate scattered light correction, while de Zeeuw et 
al.’s stays flat. Our Fe5335 profile is much steeper than the 
published ones especially outside 15 arcsec (again likely due 
to inadequate scattered light correction). Two of the index- 
combinations studied in this paper, [MgFe]’ and Tvlgb/ZTe), 
include Fe5335. The discrepancy between our Fe5335 profile 
NGC 3384 and published ones does not affect [MgFe]’ sig¬ 
nificantly. Fisher et al.’s profile (cyan curve in Fig.|HJ shows 
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Figure 3. Comparison of line strength gradients between this 
work and published values for NGC 4472. Triangles are from this 
work. Large circles are from Trager et al. (1998). Open squares 
are IFU data from Peletier et al. (1999) and crosses are long-slit 
data from Vazdekis et al.(1997). At the top of each figure, error 
bars without points represent the scatter in stars from our trans¬ 
formations to the Lick system; error bars with squares indicate 
uncertainties in Peletier et al.’s data, including those involved in 
transforming to the Lick system. 
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Figure 4. Comparison of line strength gradients between this 
work and published values for NGCs 7332, 3384, and 7457. Tri¬ 
angles are from this work. Asterisks are azimuthally averaged 
IFU data from Sil’chenko (1999), Sil’chenko et al. (2003), and 
Sil’chenko et al. (2002) for NGCs 7332, 3384, and 7457 respec¬ 
tively. Squares are SAURON data extracted along the major-axis 
from Falcon-Barroso et al. (2004), de Zeeuw et al. (2002), and 
Sil’chenko et al. (2002). For NGC 3384, crosses are long-slit data 
from Fisher et al. (1996). 


good agreement with ours. For Mgb/(Fe) however, we ob¬ 
tain a positive gradient while Fisher et al. found no gradient 
(Fig. unji, possibly due to scattered light in our data. The 
other two objects affected by scattered light are NGCs 2787 
and 3945. NGC 2787 shows no gradient in Mgb/(f’e). NGC 
3945 has an asymmetric positive gradient. This could not 
be due entirely to scattered light since the scattered profile 
was symmetric. 

The third column of Fig. IH compares profi l es for 
NGC 7457 with IFU data froni lsil’chenko et alJ (|2003) 
and archival SAURON data which were also presented in 
Silchenko et al.’s paper. Silchenko et al.’s values are system¬ 
atically smaller in the central regions than those obtained by 
others except for the case of Fe5335, where they are in agree¬ 
ment with ours. Our values have larger scatter but are other¬ 
wise in agreement with the SAURON values. The agreement 
among studies is worst for H/3 although this galaxy has lit¬ 
tle or no Balmer emission. The 2D SAURON map shows a 
negative gradient along the major axis but not along the 
minor axis while Sil’chenko et al.’s 2D map shows no gra¬ 
dients whatsoever; our values are in better agreement with 
SAURON’s than Sil’chenko et al.’s. 


3.7 Absorption Line Indices and SSP Models 

We have measured all 25 Lick indices in our galaxies as a 
function of galactocentric radius. In this paper, we concen¬ 
trate on a subset of these indices which are sensitive to age, 
metallicity, and a/Fe. The most age-sensitive indices are the 
Balmer indices H/3, Hy^, H 7 F, H(5a, and Of these, H/3 
suffers most from line-filling by emission while the Hd in¬ 
dices are the least affected. On the other hand, H/3 offers 
the most orthogonality with respect to metallicity-sensitive 
indices. Using a combination of the Balmer indices, we can 
obtain more reliable age estimates than with just one in¬ 
dex. For metal lines, we compute the indices Mgb/(Fe) and 
[MgFe]’ as discussed in TMB; the former is directly related 
to a/Fe, and [MgFe]’ traces metallicity without any sensi¬ 
tivity to a/Fe. Individually, [MgFe]’ and the Balmer indices 
are degenerate in age and metallicity but together they can 
break the degeneracy since each index has a different age- 
met allicity dependence. 

The integrated-light spectrum of an object is a linear 
combination of SSPs. Some objects, such as globular clus¬ 
ters, are well represented by a single SSP while galaxies are 
generally not. Still, one can characterize the SPs of a galaxy 
by an “equivalent SSP”. Since the integrated light from a 
galaxy is weighted by luminosity, its SSP age and metallic¬ 
ity is likely to be different from the mass-weighted age and 
metallicity of its stars. SSP values are useful parameteri- 
zations of the SPs but cannot be interpreted as true ages 
and metallicities, since galaxies most likely contain a range 
of both. To avoid over-interpreting our data, we focus on 
the line strengths, mentioning SSP values only to illustrate 
dramatic differences between objects or regions within an 
object (i.e. 2 vs 10 Gyr as opposed to 8 versus 12 Gyr). 
Different line strengths in different objects (or within an ob¬ 
ject) imply different SPs; the SSP models allow us to infer 
the underlying source of the differences. 
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3.8 Bulge-to-disk decomposition 

One-dimensional and two-dimensional bulge-to-dlsk decom¬ 
position (B/D) was performed on our images. The disk and 
bulge were simultaneously fit with an exponential and an 
exp[r^^"'] profile, respectively. Initial fits were made using 
fixed n of 1, 2, 3, and 4; the best fit of these was used as 
a starting guess in a final fit where n was allowed to vary 
as a free parameter. The results from the ID decomposition 
were used as starting guesses for the 2D. In the 2D fits, both 
the bulge and disk components were allowed to be elliptical. 
This resulted in bars being fit as bulges in some galaxies (IC 
302 and NGCs 266, 765, 3681, 3883, 5375) but not others (IC 
267 and NGCs 2487 and 5020). A two-component bulge/disk 
model, such as ours, also overestimates the bulge component 
of our three SOs with luminous inner disks (NGC s 2787 , 
3384, and 3945) as previously found bv lErwin et al.l i200,3ll . 
We pay close attention to these effects when studying the 
line strength profiles in the bulge- and disk-dominated re¬ 
gions (Section 4.2). It was found through visual inspection 
that the best-fitting bulge and disk components of the other 
galaxies were reasonable. Fig. shows the light profiles of 
our galaxies along with the best-fit bulge and disk compo¬ 
nents. 

The bulge-to-disk decomposition used here primarily 
serves to determine the relative contribution of bulge and 
disk light as a function of radius; while B/D decomposition 
is notoriously difficult, especially when allowing for a Ser- 
sic bulge, this ratio is determined more robustly than the 
derived values of bulge effective radius and Sersic index. 

Throughout this paper, the term “disk” refers to the 
exponential outer region of the luminosity profile and the 
term “bulge” refers to the excess light in the inner regions 
that was fit with a Sersic profile. It is possible that these 
photometrically-determined components do not correspond 
to a cold and a hot component, a fiat and a spheroidal com- 
pone nt, or a young and a n old component. Through simula¬ 
tions, |Abadret_^ (I 2 OO 3 II found that while the stars of a sim¬ 
ulated galaxy were well fit by a Sersic -|- exponential profile, 
the hot and cold components were not well fit individually 
by a Sersic and an exponential profile respectively. Instead 
both were Sersic in the inner regions and exponential further 
out. Using our kinematic information, we address this issue 
in Paper If. The structural properties and how they relate 
to the SPs will also be presented in that paper. 

4 RESULTS 

4.1 Central Line Strengths 

Fig. El shows the measured central values of H/3 and 
[MgFe]’ for o ur sample as well as publi shed values for th e 
Milky Way ( Rizi^^^Mjj200^ M31 llPuzia et alJ 1200^ . 
oth er bulges Proctor_lL.SarLSonj||2()()^ and ellip tical galax¬ 
ies (iTrager e^nr^9^ r ^rocto^^^msomll200^ . Our cen¬ 
tral indices were measured on spectral extractions binned to 
match Trager et al.’s 4 arcsec aperture and Proctor et al.’s 
3.6 arcsec aperture as closely as possible (3.3, 4.2, and 3.8 
arcsec for DIS I, If, and Iff respectively). Symbols are as in 
Fig.m Larger symbols denote larger central velocity disper¬ 
sions. If a point has an accompanying vertical line segment, 
its location was determined using an average of the other 


four Balmer indices instead of H/3 to compute the SSP age. 
The other end of the line segment shows the H/3 value. If the 
difference between these is large, this is due most likely to 
errors in H/3 emission correction since the galaxies whose H/3 
values lie outside the model grid have strong H/3 emission. 
Two galaxies have H/3 values that put them outside the plot 
range. These are NGCs 2787 and 5719 with H/3 values of 
-0.506 and 0.607 respectively. 

Bulges can be divided into three groups according to 
which region of Fig. El they populate: the old metal-rich 
(OMR, lower right) region , the young metal-rich (YMR, 
upper right) region with ages less than 3 Gyr and super¬ 
solar metallicity, or the metal-poor (MP, left) region with 
sub-solar metallicity. This classification scheme is analogous 
to that of Prugniel et al. with our MP, OMR, and YMR 
classes corresponding to their A, B, and C classes respec¬ 
tively. Membership in a region is closely related to bulge 
color. Most of the red bulges populate the OMR region while 
all but one of the bulges in the MP region are blue. The ex¬ 
ception is IC 267 (the red filled square with [MgFe]’=1.65). 
This bulge has strong Balmer emission which suggests that 
its red color is due to dust from recent star formation as op¬ 
posed to an old SP. Hereafter, when we discuss blue bulges, 
we include IC 267 since it is similar, in most respects, to 
blue bulges and since its B-K color puts it right on the bor¬ 
derline between blue and red classes anyway. Color does not 
uniquely represent the SPs of blue bulges; some bulges are 
blue because they are metal-poor while others are blue be¬ 
cause they are young. To distinguish the two classes of blue 
bulges in subsequent plots, the MP bulges are marked with 
an additional blue dot. 

Central line strengths are also sensitive to Hubble type. 
All the early-type (SO-Sab) bulges have large central metal- 
licities. The blue early-types are in the YMR region while the 
red early-types are in the OMR region. Bulges of late-type 
spirals are more heterogeneous in their central line strengths 
than those of early-types. The MP region is populated ex¬ 
clusively by late-types but late-types are also found in the 
other two regions. Trager et al.’s ellipticals mostly populate 
the OMR region with a couple of them lying in the YMR 
region. Lower-luminosity ellipt icals span a similar range of 
SSP values as our blue bulges iCaldwell et al.l^200,3^ : unfor¬ 
tunately, we cannot include Caldwell et al.’s data in Fig. El 
because they did not measure the same indices. The bulges 
of the MW and M31 lie in the MP region. 

Much of the variation in central line strengths is due to 
correlations between the line strengths and the global kine¬ 
matics. The different types of bulges (MP, YMR, and OMR) 
and ellipticals form a continuous and overlapping sequence 
on a plot of [MgFe]’ versus central velocity dispersion (Fig. 
Q top-left). In both bulges and ellipticals, [MgFe]’ is cor¬ 
related with (To at the low-(Jo end. As (tq increases beyond 
log (Jo >2.2, [MgFe]’ remains constant. Bulges show larger 
scatter than ellipticals in the [MgFe]’-(Jo relation. 

The apparent dichotomy seen in Fig. El between 
the MP bulges and the other objects (namely that the 
MP bulges have smaller central [MgFe]’) is naturally 
explained by the wide range of (Jq and Wmax values 
spanned by the different types of objects. The MP bulges 
reside in considerably shallower potential wells than the 
bulges and ellipticals which populate the OMR and YMR 
regions. It is important to note here that the predomi- 
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Figure 5. Results of 2D bulge-to-disk decomposition. The light profile of each galaxy is shown as a solid line, the best-fit bulge as a 
dotted line, and the best-fit disk as a dashed line. The position angles are the same as those used for the spectroscopy (Table 3). 


nant view that ellipticals are metal-rich l|Tra^er et al|2000^: tational velocity (Fig. |3 top-right), as previously found by 


nificant outliers in the [MgFe]’-Vmax, having smaller central 
values of [MgFe]’ than their red counterparts. 


lEisCTistein^jfogg^Fukngitm Nj^amura ^Bernurdi. Finkbeiner. SchBhgrfnBl'ialkaLailih.eCMnndildgfeBMi.lf&tlfa^lv&z2i2Iaagifain'T;Ovedav. Munn. I 

1200 ^ 106111001 X 6 ^al.ll200^lThcmas et ahll^nn^ applies only 

to high-cTo ellipticals. All of Trager et al.’s ellipticals have 
(To > 100 km s“^. The three ellipticals from Trager et al’ 
with the smallest velocity dispersions (100 < ao < 105 km 
s“^) include the two with the smallest SSP metallicity and 
the one with the smallest SSP age. Caldwell et al. (2003) 
found that ellipticals with ao < 100 km s“^ span a similar 
range of SSP ages and metallicities as our blue bulges. 


Balmer indices are anti-correlated with ao and weakly 
anti-correlated with Vmax (middle panel of Fig.|7|l. Residu¬ 
als in the H/3 —ctq and YiJ—Vmax relations are correlated with 
color such that at a given ao and Ymax, blue bulges (both 
MP and YMR) have larger H/3 values than red bulges. 


[MgFe]’ is also correlated with the maximum disk ro- 


While individual indices are degenerate in age and 
metallicity, the different symbol sizes in Fig. allow us 
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Figure 6. H/3 versus [MgFe]’ in the central regions of bulges and ellipticals. Magenta crosses are galaxies from Trager et al. (1998) that 
were classified as elliptical and had H/3 uncertainty less than 0.2 and velocity dispersion uncertainty less than or equal to 10 km s“^. 
Green symbols are bulges and ellipticals from Proctor & Sansom (2002). The Milky Way (Puzia et al. 2002) and M31 (Puzia et al. 2005) 
are shown as ’+’ symbols. Symbols for our bulges and those of Proctor et al. are as in Fig.0 Larger symbol sizes were used for objects 
with larger central velocity dispersion. Bulges for which we do not have color information are in black. The metal-poor (MP) bulges are 
marked with an additional blue dot. Where spectra exist along both minor and major axes, the average of the two values was used. If a 
point has an accompanying vertical line segment, its location was determined using an average of the other four Balmer indices, instead 
of H/3, to compute the SSP age. The other end of the line segment shows the H/3 value. TMB’s models are superimposed on the plot. 


to determine how metallicity and age vary with (tq- Mean 
SSP metallicity and SSP age are larger in larger -(jq galax¬ 
ies but at fixed (tq the scatter in SSP metallicity and SSP 
age are large. At fixed (tq, age and metallicity are known 
to be anti-correlated i n ellipticals llTrager et al.l l2nnnbt 
IProctor fc Sansomll20n2|) . The ellipticals from Trager et al. 


with large SSP age have smaller SSP metallicity than those 
with small SSP age. Such an anticorrelation appears to ex¬ 
ist for red bulges but with considerably larger scatter. In¬ 
terestingly, the small-(To blue bulges also appear to have an 
age-metallicity anticorrelation. 

In bulges and ellipticals, the ajFe ratio as indicated by 




























12 Bhasker K. Moorthy & Jon A. Holtzman 



1.8 2 2.2 2.4 2.6 1.8 2 2.2 2.4 2.6 


cr 


0 


V 


max 


Figure 7. Central line strengths versus central velocity dispersion and maximum disk rotational velocity. Symbols as are in Fig. 
Arrows represent upper limits for galaxies whose velocity dispersions are close to or below our resolution limit. IC 302 and NGCs 765, 
2487, 2916, and 6246A are not shown on the right panel since their inclinations are too low to measure rotation. In the bottom panel, 
the region within the inner horizontal lines corresponds to models with solar a/Fe for metallicities -1.35 ^ [Z/H] ^ 0.35 and ages from 
8 to 15 Gyr (from Fig. 4 of TMB). The region within the outer horizontal lines represents models with the same metallicities and ages 
from 3 to 15 Gyr. 


Mgb/(Fe), is correlated with ctq and Ymax (bottom panel 
of Fig. 0. In these two plots, the region within the inner 
horizontal lines corresponds to models with solar a/Fe for 
metallicities -1.35 ^ [Z/H] ^ 0.35 and ages from 8 to 15 Gyr; 
the region within the outer lines represents models with the 
same metallicities and ages from 3 to 15 Gyr (from Fig. 4 of 
TMB). Red bulges and ellipticals show good overlap in the 
Mgb/(Fe)-o'o diagram. With a few exceptions, blue bulges 


are consistent with having solar a/Fe. Consequently, most 
blue bulges have smaller Mgb/(Fe) ratios than red bulges 
or ellipticals at a given value of (tq or Ymax- One of the 
blue bulges, NGC 6246A, has super-solar alFe, large age, 
and small metallicity ([Z/H]—0.8) like MW halo stars. The 
other two blue bulges with super-solar ajFe have super¬ 
solar metallicity like the majority of red bulges. 

There are hints that barred galaxies follow different 
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Figure 8. [MgFe]’ profiles in our large-bulge galaxies. Open squares and filled triangles show points that were and were not corrected for 
[NI] 5199 emission respectively. For NGC 3384, the cyan curve shows the profile obtained by Fisher et al. (1996). At the top-left of each 
plot are symbols denoting the galaxy type (as in Fig.|3 as well as the NGC or IC identifiers. The solid and dotted vertical lines indicate 
where the ratio of bulge to disk light is two and half respectively, as determined from the bulge-to-disk decomposition. The solid and 
dotted red arrows indicate the location of the bulge effective radius and disk scale length respectively. Results from linear least-squares 
fits performed separately in the bulge- and disk-dominated regions are shown in red. Blue lines are color profiles (B-K-3) from PB and 
DJ. B-K was not available for IC 302; V-H-3 is shown instead. 

index-fjo and index-V^aa: relations than unbarred galaxies. 

At fixed CTO and Vmaa, barred galaxies appear to have larger 
central values of [MgFe] ’ than unbarred galaxies (or galaxies 
with elliptical-shaped bulges) of the same ctq or Vmax, with 
IC 267 being the most notable exception. B/p bulges gener¬ 
ally lie between and exhibit larger scatter than barred and 
unbarred/elliptical bulges in the [MgFe]’-CTo and [MgFe]’- 


Vmaa: diagrams. The central regions of b/p bulges could be 
contaminated by the foreground disk resulting in smaller 
values of [MgFe]’ than those of the low-inclination barred 
galaxies. Barred galaxies appear to have smaller H/3 than 
unbarred galaxies at fixed ctq but not at fixed Ymax- For 
2.2 > CTO > 2.35, barred galaxies appear to have larger val¬ 
ues of Mgb/(Fe) but the opposite is true for 2 > ctq > 2.2. 
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Figure 9. [MgFe]’ profiles in our small-bulge galaxies. Symbols 
are as in Fig. 

No striking difference is seen between barred and unbarred 
galaxies in the Mgb/{Fe)-Vmaa: relation. These results need 
to be verified using a large sample of barred and unbarred 
galaxies as a function of no and Vmax- 

The central regions of our three barred SOs with disk¬ 
like structure and kinematics (filled black circles) have 
super-solar SSP metallicities and a/Fe ratios and two of 
the three have large SSP ages. If all disks formed stars on 
long timescales (several Gyr), we would expect them to have 
small a/Fe ratios, like the MW disk at the solar neighbor¬ 
hood, since ISM enrichment would eventually be dominated 
by SN la. The “luminous inner disk s” of these SOs have 
not had such a star-formation history. IPeletier et alJ lll99flll 
found that the Sombrero galaxy is also dominated by a fast 
rotating disk whose [Mg/Fe] is similar not to other disks but 
to ellipticals of similar mass as the Sombrero. 

4.2 Line Strength Gradients 

Figs. ISI9I show spatially resolved [MgFe]’ profiles in our 
galaxies. They are split into two figures according to whether 
they have large or small bulges. Open squares and filled tri¬ 
angles show points that were and were not corrected for [NI] 
5199 emission respectively; the correction was seldom per¬ 
formed. The solid and dotted vertical lines indicate where 
the ratio of bulge to disk light is two and one half respec¬ 
tively, as determined from the bulge-to-disk decomposition. 
In the central regions, indices were measured on approxi¬ 
mately 1 arcsec bins. 

The large-bulge galaxies have five or more measure¬ 
ments within the solid vertical lines while the small bulge 
galaxies have three or fewer. We performed linear least- 
squares fits to the [MgFe]’ profiles separately in the bulge- 
and disk-dominated regions, selecting the points to include 
based on the bulge-to-disk decomposition but making excep¬ 
tions where they seemed appropriate (such as when the bar 
was fit as a bulge; see Section 3.8). We usually avoided the 


Figure 10. The [MgFe]’ gradient within the bulge-dominated 
region versus the central value of [MgFe]’. Symbols are as in Fig. 

m 



Figure 11. The central value of [MgFe]’ versus the value at one 
disk scale length computed using the results from a least-squares 
fits to the [MgFe]’ profiles. The solid line shows where objects 
would lie if their bulges and disks were identical in [MgFe] ’. Sym¬ 
bols are as in Fig.|S| 

transition region from bulge to disk dominance. The best-fit 
lines are shown in red. The large-bulge galaxies have neg¬ 
ative [MgFe]’ gradients (weaker [MgFe]’ with increasing ra¬ 
dius) in the bulge-dominated region. In all but two of these, 
[MgFe] ’ decreases steadily from the galaxy center to the solid 
vertical lines, beyond which the slope of the profile changes. 
The exceptions are NGCs 3681 and 5707, where [MgFe]’ is 
constant within the solid vertical lines but decreases in the 
bulge-disk transition. In low-inclination galaxies and along 
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Figure 13. H<5a profiles in our large-bulge galaxies. Symbols are as in Fig. El 


the major axes of inclined galaxies, [MgFe]’ is usually larger 
just outside the solid vertical lines than inside it (NGC 5987 
is a good example). The fact that we can identify distinct 
bulge and disk components in the [MgFe]’ profiles suggests 
that disk contamination is not significant within the solid 
vertical lines. 

Another test for disk contamination is how the major 
and minor axis profiles vary as a function of inclination. Low- 
inclination galaxies should have identical profiles if there are 
no azimuthal differences in line strengths. The major and 
minor axis profiles of NGC 2775, the only low-inclination 
galaxy for which both were obtained, are indeed identical. 


The minor axis profiles of the edge-on large-bulge galaxies 
(NGCs 5422, 5689, 5746, 5793) continue decreasing beyond 
the distance at which the m ajor axis profiles flat ten off, a re¬ 
sult previously obtained bv iFisher et alJ (Il99fil for edge-on 
SOs. This implies that the outer bulge has lower metallicity 
than the inner disk. At intermediate inclinations, one ex¬ 
pects more disk contamination on one side of the minor axis 
(the dusty side) than the major axis for the same solid an¬ 
gle. In thethree intermediate inclination galaxies for which 
we have major and minor axis spectra (NGCs 3544, 5326, 
and 5389), this effect is clearly seen; the profile is asymmetric 
outside the solid vertical lines with the profile turning over 
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Figure 12. [MgFe]’ gradient versus B-K color gradient, where a 
gradient is defined as the difference between the value at the cen¬ 
ter of the galaxy and the value at the disk scale length. Symbols 
are as in Fig. 


at a smaller galactocentric distance on the dusty (left) side 
than on the dust-free side. Within the solid lines, there is 
good agreement in gradient slopes between major and minor 
axes, indicating that disk contamination is not significant. 

The small-bulge galaxies generally do not have nega¬ 
tive [MgFe]’ gradients within the bulge-dominated region. 
In four of them (IC 267 and NGCs 3831, 5577, and 6246A), 
there is the hint of a positive gradient while the rest (NGCs 
2916, 5362, and 6368) are consistent with having little or no 
gradient. Since the profiles are always different inside and 
outside the solid lines, they cannot be explained by disk 
contamination. 

The slopes of the [MgFe]’ gradients within the bulge 
are shown in Fig. uni as a function of the central [MgFe]’ 
values. Galaxies with large central values have correspond¬ 
ingly large negative gradients while the three galaxies with 
the smallest central values, namely IC 267 and NGCs 5577 
and 6246A, have positive gradients. Since central [MgFe]’ 
is correlated with ctq and Ymax, it follows that the slope 
of [MgFe]’ gradients in bulges is correlated with the global 
kinematics. 

Most galaxies also have negative gradients in the disk- 
dominated region but it is shallower than that of the bulge. 
Some galaxies (e.g. NGCs 5746, 5838, and 7332] have no 
gradient in the disk. 

The [MgFe]’ value at one disk scale length (computed 
using the results of our least-squares fits) is correlated with 
the central value IFig. Hill . This indicates that the metal- 
licity of the disk is correlated with that of the bulge. This 
correlation holds for all galaxies, not just those with bars, 
blue bulges, or bulges identified as having disk-like struc¬ 
tural or kinematical properties. 

The blue lines in Figs. I8I9I are color profiles {B — K — 
2) from PB and DJ. The shapes of the [MgFe]’ and color 
profiles agree often but not always. Discrepancies occur most 
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Figure 14. profiles in our small-bulge galaxies. Symbols are 
as in Fig. El 

often in the central regions. Several galaxies (e.g. IC 267 and 
NGC 266, 2487, 2916, 3728, and 5577) show a positive color 
gradient in the central 5 arcsec; of these, only IC 267 and 
NGC 5577 have a corresponding positive [MgFe]’ gradient. 
On the other hand, the central region of NGC 6246A has a 
positive [MgFe]’ gradient but its color profile is flat, de Jong 
(1996) noted that it is not possible to identify distinct bulge 
and disk components using the color profiles. However, it is 
possible to do so using the [MgFe]’ profiles. As mentioned 
earlier, the slopes of the [MgFe]’ profiles are almost always 
distinct inside and outside the bulge-dominated region, with 
the former usually having a steeper negative gradient. 11 out 
of the 14 DJ galaxies and 15 out of the 17 PB galaxies have 
negative [MgFe]’ gradients in the bulge-dominated region. 
The majority of the PB galaxies (12 out of 17 as opposed to 4 
out of the 14 DJ galaxies) also show a negative color gradient 
in the bulge-dominated region. The systematic difference in 
color gradients between the PB and DJ samples is likely due 
to different amounts of extinction at different inclinations 
since the PB galaxies have large inclinations while the DJ 
galaxies have small ones. 

If we compare the gradient in [MgFe]’ from the galaxy 
center to the disk scale length (computed using the fit re¬ 
sults) with the color gradient (read from the profiles), we 
find that these two quantities are not tightly correlated 
(Fig. [nil. This is most likely due to the discrepancy be¬ 
tween [MgFe]’ and color profiles in the central regions of 
some galaxies, which in turn is due most likely to the color 
prof iles being affected by dust. 

ICadotti fc dos AniosI i200 ji found a greater prevalence 
of null or positive color gradients in barred galaxies than in 
unbarred galaxies. They interpreted their result as evidence 
for gradients being erased by bar-driven mixing. We do not 
see any systematic difference between barred and unbarred 
galaxies with regard to their gradients. Also, if bars homog¬ 
enized the SPs, we might expect a smooth transition in the 
line strength profiles from the bulge to the bar. However, 
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Figure 15. Gradients in and [MgFe]’. Red arrows are drawn from the galaxy center to the edge of the bulge-dominated region on 
either side of the galaxy. Blue arrows are drawn from the edge of the bulge-dominated region to the disk scale length. Symbols are as 
in Fig. TMB models with solar a/Fe, age=l, 2, 3, 6, and 15 Gyr, and [Z/H]=-1.35, -0.33, 0, 0.35, and 0.67 are overlaid. Unlike H/3, 
the higher order Balmer indices are not independent of a/Ye. a-enhanced models are parallel to the solar models, lying above and to the 
right. The cyan curves show models with age=l Gyr and [Z/H]=0.67 for [a/Fe]=0.3. 


the outer bulges of barred galaxies have lower [MgFe]’ than 
the inner bar the same way the outer bulges of unbarred 
galaxies have have lower [MgFe]’ than the inner disk. 

4 ..2.1 Separating Age and Metallicity Effects 

Fig.cni shows gradients in the H5 a index. The H(5 profiles 
show less scatter than the lower-order Balmer indices since 
they are less affected by emission. Squares and filled tri¬ 


angles show points that were and were not corrected for 
emission, respectively. Least-squares fits were performed on 
the Balmer indices exactly in the same manner as for the 
[MgFe]’ profiles. The fit results for HJa are shown in red 
in Fig. uni Gradient slopes computed on individual indices 
were combined to disentangle the effects of age and metal¬ 
licity. This is shown in Fig. FTTl for the [MgFe] ’-Htf a index 
combination. A red arrow is drawn from the galaxy center 
to the edge of the bulge-dominated region (on either side of 
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Figure 16. Mgb/(F’e) profiles in our small-bulge galaxies. For NGC 
(1996). Symbols are as in Fig.|^ The horizontal lines are as in|7| 

the galaxy) and a blue arrow is drawn from there to the disk 
scale length. 

The majority of galaxies (at least 29 out of 38) have 
negative metallicity gradients in the bulge-dominated re¬ 
gion. NGCs 3681, 3831, 5362, 5707, and 7311 show little 
or no metallicity gradient in the bulge. Three of the five MP 
bulges (IC 267 and NGC 5577 and 6246A) have positive 
metallicity gradients. The remaining two galaxies (NGCs 
5746 and 5793, both edge-on with b/p bulges) show inter¬ 
nal discrepancies in the fit results. Except for its minor axis 
H(5 a profile, NGC 5746 is consistent with having little or no 


3384, the cyan curve shows the profile obtained by Fisher et al. 


metallicity gradient. The minor axis profiles of NGC 5793 
consistently a negative metallicity gradient but the major 
axis profiles show none. 

The majority of galaxies are consistent with having lit¬ 
tle or no age gradient within the bulge. At least ten galaxies 
(IC 302 and NGCs 266, 765, 2487, 2599, 2916, 3883, 5020, 
5375, and 5838) have a positive age gradient (larger age with 
increasing radius). Of these, seven are barred and one has a 
b/p bulge. At least five galaxies (NGCs 3728, 5577, 6246A, 
7311, and the major axis of NGC 5793) have a negative age 
gradient in the bulge. 
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Figure 17. Mgb/(Fe) profiles in our small-bulge galaxies. Sym¬ 
bols are as in Fig.|^ The horizontal lines are as in|2] 
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Figure 18. Bulges on the BPT (Baldwin et al. 1981) diagram 
in which the emission line flux ratio [0111]5007/ H fS is plotted 
against the ratio [NII]G583/Ha. Symbols are as in Fig. The 
dashed curve shows the demarkation between starburst galaxies 
and AGN as defined by Kauffmann et al. (2003). We adopt the 
commonly used definition of LINERs having [OIII] 5007/H/3 < 3 
and[NII] 6583/H« > 0.6 and Seyferts having [OIII] 5007/H/3 > 3 
and[NII] 6583/Ha > 0.6. 


The Balmer indices show large scatter in the disk- 
dominated regions, making it difficult to determine whether 
metallicity and age gradients are present. Plots of gradients 
in the five Balmer indices versus [MgFe]’ (such as Fig. I15II 
do not show consistent results for approximately half the 
galaxies. Of the other half, most have a negative age gra¬ 
dient and little or no metallicity gradient. Most disks are 


solar or sub-solar in metallicity but NGCs 5746 and 5838 
are super-solar well into the disk-dominated region. 


4-2.2 Abundance Ratio Gradients 

Most galaxies either have a positive gradient or no gradient 
in Mgb/(F’e) within the bulge-dominated region (Figs. 1161 
and IT7t . The disk-dominated regions generally have solar 
a/Fe. Since the red-bulge galaxies have super-solar a/Fe in 
the center, they have a negative gradient in the bulge-disk 
transition. The blue-bulge galaxies have solar a/Fe in the 
center. These either have uniformly solar a/Fe or a positive 
gradient in the bulge and a negative gradient in the bulge- 
disk transition. Recall that what is marked as the bulge- 
dominated region in some galaxies (e.g. NGCs 266 and 5375) 
is actually a bar and that the true bulge-dominated region 
is smaller. The elliptical galaxy, NGC 4472, is uniformly 
super-solar in a/Fe. 

There are a few galaxies that have super-solar a/Fe 
in the disk-dominated region. NGCs 266, 5707, and 5746 
are nearly uniformly super-solar. The disk of NGC 5838 is 
super-solar but less enhanced than its bulge. 


4.3 Emission Lines in Bulges 

Fig ■ llOl shows profiles of Ha and [Nil] 6583 emission strength 
in our galaxies. We detect emission in the central regions of 
all our galaxies except the SOs NGCs 3384 and 7457. The 
locations of our galaxies on the BPT iBaldwin et al.lli98lll 
diagram of emission line ratios is shown in Fig. ll8lfor objects 
with central emission-line EW smaller than -0.5 A (the nega¬ 
tive sign denotes emission) in Ha, [Nil] 6583, H/3, and [OIII] 
5007 . The dashed curve shows the dem arkation between 
starbu rst galaxies and AGN as defined bv iKauffmann et alJ 
i2003h . The majority of our emission-line galaxies are AGN. 
If the emission in these galaxies is due entirely to the AGN, 
we would expect it to be restricted to the center of the 
galaxy. However, in the three of these (NGCs 2599, 5719 
and 7537), it is not centrally peaked. In the other five (NGCs 
266, 2916, 3831, 5793, and 6368), the emission peaks at the 
center and decreases steadily out to the edge of the bulge- 
dominated region, beyond which it rises again. Therefore, 
all or most of the AGN also have active star formation in 
the bulge-dominated region. 

AGN have previously been found to ha ve a larger frac¬ 
tion of you ng stars than quiescent galaxies jRaimann et alJ 
I 2 OOII l2003n . In agreement with these results, we find that 
most of the AGN have small SSP ages (< 4 Gyr). The 
only one with a large SSP age (i5 Gyr) is NGC 3831. This 
could be due to errors from emission correction or from 
the young componen t not dominating the total luminosity. 
IPrueniel et al.l [1200 found that bulges with emission were 
small and metal-poor. The star-forming region of our BPT 
diagram is populated by four blue bulges. They have similar 
stellar populations as Prugniel et al.’s emission-line galaxies 
except that one of them (NGC 6246A) has a large SSP age, 
again possibly due to errors in emission correction. 

We see a wide range of behaviors in the emission-line 
profiles. In some galaxies, such as NGCs 3681 and 5362, 
there is strong emission in the disk-dominated region but 
little or no emission in the bulge-dominated region as would 
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Figure 19. Profiles of emission line strengths in bulges. Triangles are Ho. Squares are [Nil] 6583. Symbols are as in Fig. El 


be expected if disks continue to form stars while bulges 
do not. The only galaxies with little or no emission in the 
disk-dominated region are SOs. In other cases, there is emis¬ 
sion throughout the galaxy but it is weaker in the bulge- 
dominated region (e.g. NGCs 1642, 2916, 5020, and 7537). 
This is consistent with a quiescent bulge and a star-forming 
disk coexisting in the central regions with the ratio of bulge 
to disk dominance decreasing with radius. Alternatively, the 
bulge and disk could both be forming stars but the disk more 
actively so. Finally, in some cases (e.g. IC 267 and NGCs 266 
and 5793), the emission lines are strongest at the center. 


5 THE FORMATION OF BULGES 

As mentioned in the introduction, present-day ACDM cos¬ 
mology argues against the monolithic collapse scenario as 
does observational evidence for the recent and continuing 
mass assembly of ellipticals. Of the main proposed forma¬ 
tion scenarios, that leaves mergers and secular evolution as 
possibilities for bulges. 

However, the collapse model continues to receive much 
attention under the claim that it better reproduces the ob¬ 
served line strength profiles of ellipticals. We investigate 
whether or not this is true for bulges. Gradients in the 
index Mg 2 have been computed in galaxies formed in col- 
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Figure 20. Comparison of Mg 2 profiles between this work and numerical simulations. Solid lines are merger models by Bekki & Shioya 
(1999) final age 13.1 Gyr and initial disk masses of IO^'^Mq (bottom curve) and IO^^Mq. Dotted lines are collapse models by Angeletti 
&: Giannone with final age 13 Gyr (top curve) and 2 Gyr. Symbols are as in Fig. 


lapse and merger simulations, allowing for direct compar¬ 
isons with our data (Figs. I2(H . Points are our data. Solid 
lines are two remnan ts from major disk-disk mergers by 
iRekki fc Shioval ill 99^1 with initial disk masses of 10 ^°Mq 
( bottom curve in Mg 2 ; top curve in H/3) and W^^Mq. The 
remnants are 13.1 Gyr old. Dotted line s are two collapse 
models by lAngeletti fc Oiannonel ll2^^.3^ with final ages of 
13 Gyr (top curve) and 2 Gyr. We focus on the gradient 
slopes, assuming that changes in mass and formation epoch 
mainly shift the models up or down. The collapse models 


predict steeper Mg 2 profiles than the merger models. Within 
the bulge-dominated region, some of our profiles agree with 
one of the merger models (e.g. NGCs 765 and 7311) while 
others agree with one of the collapse models (e.g. 3728 and 
the minor axes of NGCs 5326 and 5422). There are also 
cases, mostly among blue bulges, where the observed pro¬ 
file is flatter than the merger models (e.g. IC 302). However, 
the majority of galaxies fall between the collapse and merger 
models. 

Bekki & Shioya also computed H/3 profiles in their mod- 
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merger models by Bekki & Shioya (1999) with initial disk masses of 


Figure 21. Comparison of H/3 profiles between this work and 

Mq (bottom curve) and 10^^ Mq. Symbols are as in Fig.lHl 

els (Fig. I2H . Nearly all our galaxies have flat H/3 profiles 
within the bulge-dominated region as predicted by the mod¬ 
els. The profiles of the oldest bulges agree with the models 
in their zero-points as well, while younger bulges lie above 
the models. 

Through chemical evolution modeling, Thomas et al. 
(1999) studied a/Fe ratios in ellipticals that formed through 
a fast (1 Gyr) collapse of star-forming clumps and through 
mergers of MW-like spirals. The main difference between the 
two models was that the merging spirals had several Gyr of 
Fe enrichment while the gas involved in the collapse was 


not pre-enriched in Fe. The large central a/Fe ratios found 
in massive ellipticals were reproduced in the collapse model. 
The merger model does not produce super-solar a/Fe assum¬ 
ing a Salpeter IMF unless the merger happened early, before 
the progenitors acquired much Fe. Metallicity and a/Fe are 
anticorrelated in the collapse model; since ellipticals have 
negative gradients in metallicity, the model predicts that 
they should have a positive gradient in a/Fe. The merger 
model produces solar a/Fe in the outer regions. Therefore, 
uniformly solar a/Fe is consistent with the merger model 
while positive gradients are consistent with the collapse 
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model. Pipino et al. (2005) also find positive a/Fe gradi¬ 
ents in a model where ellipticals are formed through the in¬ 
fall of gaseous lumps. Since our blue bulges have uniformly 
solar a/Fe, they are consistent with the predictions of the 
merger model. Several of the red bulges have positive gra¬ 
dients as predicted by the collapse model. Most of the re¬ 
maining bulges have uniformly super-solar a/Fe, which is 
difficult to reproduce in any of the models. 

In summary, both collapse and merger models have lim¬ 
ited success in reproducing the line strength proHles of in¬ 
dividual galaxies but neither explains the full range of be¬ 
haviors seen in the data. It is important to note that hier¬ 
archical models are only beginning to make robust predic¬ 
tions for SPs, successfully reproducing properties tradition¬ 
ally thought to favor the collapse model, such as the mass- 
metallicity relation. As advancements continue to be made 
in incorporating gas dynamics, star formation, and chemi¬ 
cal evolution in cosmologically-motivated merger models, it 
will be interesting to see if the line strength profiles will be 
reproduced as well. 

5.1 Mergers 

In a recent paper, iFaber et al.l (1200, argued that massive 
red ellipticals could not have formed entirely through major 
mergers of gas-rich components or through dry mergers but 
through a combination of the two. Ellipticals of the same 
mass and color could have formed in different ways: through 
early gas-rich mergers of low-mass objects followed by dry 
mergers or through recent gas-rich mergers of more massive 
objects. Objects that arrived on the red sequence early-on 
and have been gaining mass through dry mergers will have 
larger SSP ages than those that have arrived on the red 
sequence near their present mass as the result of recent gas- 
rich mergers. The former will also have smaller metallicities 
since their last gas-rich mergers were of lower mass progen¬ 
itors with correspondingly lower metallic ities according to 
the gas-phase mass-metalli city relation iKobulnickv et ^ 
1200.4 iTremonti ~ IDB), The predicted anticorrelation 
between age and metallicity at fixed no is seen in ellipticals. 
If such an anticorrelation exists for bulges, it is not nearly as 
tight as that of ellipticals. This suggests that an additional 
formation mechanism might be required to explain the SPs 
of bulges. 

Recent semi-analytic models incorporating the Mille¬ 
nium Simulation of cosmic structure growth find a cor¬ 
relation between stellar metallicity and stellar mass, with 
the most massive gal axies having roughly solar metallicity 
JPe Lucia et al.ll2005ll. This is qualitatively consistent with 
the observed [MgFe]’-(To and [MgFe]’-Vmaa: relations. Mas¬ 
sive ellipticals and bulges have super-solar central metallic¬ 
ity which is in apparent contradiction with de Lucia et al.’s 
results. However, they also have negative metallicity gradi¬ 
ents. The arrows in Fig. EHl which extend to approximately 
the bulge effective radius, fall around solar metallicity in 
the massive red bulges. The metallicity at the effective ra¬ 
dius is more representative, than the central value, of the 
mean metallicity. Therefore, the data are not inconsistent 
with the models. 

de Lucia et al. also find in their simulations that less 
massive ellipticals had more extended star formation histo¬ 
ries than their massive counterparts. This is consistent with 


the observed a/Fe-no and a/Fe-Vmaa: correlations. These 
correlations can be produced in starbursts induced by gas- 
rich mergers. During the starburst, SN II enrich the ISM 
with a-elements. If star-formation is somehow quenched be¬ 
fore SN la contribute much Fe, the a/Fe ratio increases. Dry 
mergers would add scatter to the a/Fe-no and a/Fe-Vmai 
relations since they increase cro and Vmax without altering 
the a/Fe ratio. 

The differences in between blue and red bulges at fixed 
(To can also be explained by mergers. At fixed no, blue bulges 
have smaller SSP ages than their red counterparts which 
suggests that they have undergone gas-rich mergers more 
recently. The progenitors of the blue bulges have then had 
more time to acquire Fe. Since the progenitors have small 
a/Fe ratios, so do the remnants. This is seen in the sim¬ 
ulations by Thomas et al. (1999), who found that gas-rich 
mergers cannot produce large a/Fe ratios unless they hap¬ 
pened early in the chemical evolution of the progenitors. 

5.2 Secular Evolution 

In dissipationless secular evolution, the bulge is formed 
through the vertical and radial redistribution of disk stars. 
In this process, existing gradients can either become am¬ 
plified since the resulting [pseudo] bulge has a smaller scale 
length than the progenitor disk or erased as a consequence of 
disk heating. But if the disk has no gradient, neither should 
the bulge. This process cannot be ruled out on the basis of 
[MgFe]’ gradients since the majority of galaxies show neg¬ 
ative gradients in the disk-dominated region. However, the 
majority of red bulges have solar a/Fe in the disk-dominated 
region despite having super-solar a/Fe in the bulge. There¬ 
fore they could not have been produced through purely dis¬ 
sipationless secular evolution. If secular evolution with gas 
infall has been responsible for the formation of these ob¬ 
jects, the star-formation timescales must have been identi¬ 
cal (at fixed no) in this scenario as in merger-induced star- 
formation since red bulges and ellipticals follow the same 
a/Fe-(7o relation. Furthermore, the star formation must have 
been completed several Gyr ago since red bulges have large 
SSP ages. This goes for the three barred SOs with disk-like 
structural and kinematical properties (NGCs 2787, 3384, 
and 3945) as well. These objects are identical to ellipticals of 
comparable (tq in their stellar populations and two of them 
have among the largest central SSP ages observed. 

Note that the a/Fe ratios of the blue bulges are con¬ 
sistent with dissipationless secular evolution. Unfortunately, 
neither mergers nor secular evolution can be ruled out for 
blue bulges on the basis of their a/Fe ratios. 

Secular evolution with gas infall is supported by the 
frequency of barred galaxies with age gradients. Of the ten 
galaxies whose central regions are younger than the outer 
regions, seven are barred and one have a b/p bulge. Bar- 
driven gas infall could lead to extended star-formation in 
the central region, producing the observed age gradient. 

If bars are long-lived and the chemical imprints of sec¬ 
ular evolution are different from those of mergers, we would 
expect the bulges of barred galaxies to have different abun¬ 
dance patterns than those of unbarred galaxies. We see hints 
of such differences in index-no and index-Vmai; relations. At 
fixed (To and Ymax, barred galaxies appear to have larger 
central metallicities. 
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The metallicities of bulges and their disks are corre¬ 
lated. This is naturally explained in processes that involve 
the bulge being formed from the disk. However, this corre¬ 
lation holds for all galaxies, not just those with bars, blue 
bulges, or bulges identified as having disk-like structure and 
kinematics. Therefore either all bulges formed secularly and 
some had their bars destroyed or the other bulge/disk for¬ 
mation mechanisms also produce this correlation. 

5.2.1 Evolution of Galaxy Populations 

Small-(T bulges fall into two categories: YMR bulges with 
little or no star formation and MP bulges which are ac¬ 
tively forming stars. This suggests that the MP bulges would 
have migrated to the YMR region by the time their star- 
formation is quenched. Will this be the scenario for all metal- 
poor bulges (including that of the Milky Way) or is the ob¬ 
served anticorrelation between emission strength and metal- 
licity the result of small number statistics? Are there really 
no metal-poor bulges that do not have emission? Extending 
this type of study to large galaxy samples should shed light 
into the evolution of small-cr bulges. 

While all hve of the MP galaxies are late-types (Sb-Sc), 
three of the seven YMR galaxies are early types (SO-Sa). 
Perhaps, the mechanisms that trigger and quench the star- 
formation are also responsible for tranforming galaxies from 
late- to early-types. As the YMR bulges age, they will move 
down to the OMR region. 


6 SUMMARY 

We have studied line strengths in the bulges and inner disks 
of 38 galaxies in the local universe. Our galaxies span a 
wide range of Hubble types, central velocity dispersions, 
maximum disk rotational velocities, and inclinations. The 
low-inclination galaxies include barred and unbarred ob¬ 
jects; the edge-on galaxies include those with and without 
boxy/peanut-shaped bulges. We included several galaxies 
whose bulges were previously identified as being disk-like 
in their colors or kinematics to see if their spectral proper¬ 
ties reveal evidence for secular evolution. We use the [MgFe]’ 
index and five Balmer indices to characterize the luminosity- 
weighted metallicities and ages of the SPs and the Mgb/(Fe) 
index to characterize the a/Fe ratios. Our main results are 
the following: 

• The central regions of bulges range in SSP metallicity 
from [Z/H]=-0.8 to -1-0.7 dex and in SSP age from less than 
2 to greater than 15 Gyr. 

• The central ages and metallicities are sensitive to bulge 
color which is in turn sensitive to central velocity dispersion 
and maximum disk rotational velocity. 

• Red bulges of all Hubble types are similar to luminous 
ellipticals in their central SPs. They have large SSP ages 
and are super-solar in SSP metallicity and a/Fe. 

• Blue bulges can be separated into two classes: a metal- 
poor class that is restricted to late-types with small velocity 
dispersion and a young, metal-rich class that includes all 
Hubble types and velocity dispersions. The metal-poor blue 
bulges are actively forming stars while the metal-rich ones 
are not. Low-luminosity ellipticals exhibit a similar range of 
SSP ages and metallicities as blue bulges. 


• Luminous ellipticals and the different types of bulges 
form a continuous and overlapping sequence on diagrams of 
metallicity- and age-sensitive indices versus (tq. At fixed no, 
there is no systematic difference between bulges and ellipti¬ 
cals on these diagrams but bulges exhibit larger scatter. At 
fixed (70, age and metallicity are more tightly anticorrelated 
in ellipticals than in bulges. 

• a/Fe in red bulges is correlated with (jq and Ymax- Red 
bulges and ellipticals follow the same a/Fe-cro relation. 

• Most blue bulges (11 out of 14) are consistent with hav¬ 
ing solar a/Fe. At fixed (tq, blue bulges have smaller a/Fe 
than red bulges and ellipticals. 

• Barred galaxies appear to have larger central metallic¬ 
ities than unbarred galaxies of the same (tq and Vmax. 

• Most galaxies show a steady decrease in metallicity- 
sensitive indices with radius. The slope of the gradient is 
correlated with the central value and therefore with the 
global kinematics. The bulge- and disk-dominated regions 
are distinct in their line strength proHles, with the disks gen¬ 
erally having shallower slopes. The smallest bulges do not 
have negative line strength gradients; some of these have 
flat profiles in the central region while others have positive 
gradients. 

• There is a correlation between [MgFe]’ strength in the 
bulge and the disk. This correlation holds for all galaxies, 
not just those with bars, blue bulges, or bulges identihed as 
having disk-like structural or kinematical properties. 

• Where positive age gradients (with the central regions 
being younger) are present, they are invariably in barred 
galaxies. This suggests that bar-driven star formation has 
occurred. However, several red bulges in barred galaxies have 
large central SSP ages (although it could be younger than 
the outer regions) which means there has been no signihcant 
bar-driven star formation for several Gyr. 

• Four galaxies have super-solar a/Fe in the disk- 
dominated region. The rest are consistent with having solar 
a/Fe in the disk. 

• Objects identihed as having disk-like structural or kine¬ 
matic properties do not have noticeably different SPs than 
other bulges. They follow the same scaling relations as the 
red bulges and ellipticals and have metallicity gradients. The 
three barred SOs identihed as having bulges with disk-like 
structural and kinematic properties are also a-enhanced and 
therefore do not resemble the majority of the disks, includ¬ 
ing the MW disk at the solar neighborhood. 

• Color prohles agree frequently but not always with line 
strength prohles. Where there is a discrepancy, due likely to 
the colors being affected by dust, it is usually in the central 
regions. Consequently, color gradients (computed as the dif¬ 
ference in color between the center and a characteristic scale 
length) do not necessarily correlate with [MgFe]’ gradients, 
illustrating the value of spectroscopy. 

Overall, our results are consistent with the hypothesis 
that mergers have been the dominant mechanism responsi¬ 
ble for the formation of bulges. However, some of the ob¬ 
servations, such as the correlation between bulge and disk 
metallicity, pose signihcant challenges to the merger sce¬ 
nario. Furthermore, the possibility that barred galaxies fol¬ 
low different scaling relations than unbarred galaxies and are 
overrepresented among galaxies with age gradients supports 
the secular evolution picture. 
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Central line strengths on a statistically significant sam¬ 
ple of ellipticals and bulges of barred and unbarred spirals 
would be invaluable in determining whether more than one 
formation mechanism is required for bulges. The necessary 
data are already available in the databases of large surveys 
such as the SDSS. Spatially resolved studies on a smaller, 
representative sample, would allow for better comparisons 
between gradients in different types of galaxies. 
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Table 1. The Galaxy Sample. DJ denotes de Jong (1996), PB denotes Peletier &: Balcells (1996), and EN denotes the ENEAR survey 
(Wegner et ah 2000). Morphological types are from the NASA/IPAC Extragalactic Database. B magnitudes are from RC3 (de Vaucouleurs 
et al. 1991). Bulge and disk colors are from DJ and PB. Bulge color is defined to be the color at half the K-band bulge effective radius 
or 5 arcsec, whichever is larger. Disk color is defined to be the color at two disk scale lengths, b/a is the red major over minor axis ratio 
taken from the sources listed. The recessional velocities shown are the RC3 heliocentric velocities corrected to the Local Group according 
to Karachentsev & Makarov (1996). Where both 21cm and optical velocities were available, the 21 cm values were used. Since no RC3 
data was available for NGC 3831, the B magnitude and optical heliocentric velocity were taken from Fairall et al. (1992). The scale was 
obtained assuming Ho=70. The distance to NGC 3384 was determined by Cepheid observations (Tanvir et al. 1999). 


Galaxy 

Source 

Type 

Morph. 

ms (mag) 

(B-K)s 

{B-R)d 

b/a 

Vlg (km/s) 

Scale (kpc/arcsec) 

IC 267 

DJ 

SBb 

Bar 

13.63 

4.6 

4.24 

0.71 

3577 

0.25 

IC 302 

DJ 

SBbc 

Bar 

13.59 



0.92 

5950 

0.41 

IC 1029 

PB 

SAb 

Ell* 

13.64 

3.89 

3.77 

0.24 

2520 

0.17 

NGC 266 

DJ 

SBab 

Bar 

12.33 

4.6 

3.85 

0.94 

4908 

0.34 

NGC 765 

DJ 

SABbc 

Bar 

13.60 

4.4 

3.82 

1.00 

5117 

0.37 

NGC 1642 

DJ 

SA(rs)c 

Unb? 

13.28 

4.0 

3.26 

1.00 

4579 

0.32 

NGC 2487 

DJ 

SBb 

Bar 

13.10 

4.1 

3.53 

0.92 

4758 

0.33 

NGC 2599 

DJ 

SAa 

Unb 

13.12 

4.0 

3.8 

1.00 

4651 

0.32 

NGC 2775 

EN 

SAab 

Unb 

11.13 



0.85 

1173 

0.08 

NGC 2787 


SB0+ 

Bar 

11.77 




696 

0.06 

NGC 2916 

DJ 

SAab 

Unb 

12.42 

4.2 

3.59 

0.74 

3618 

0.25 

NGC 3384 


SBO- 

Bar 

10.63 




735 

0.06 

NGC 3544 

EN 

SABa 

Unb? 

12.99 



0.30 

3354 

0.23 

NGC 3681 

DJ 

SAB(r)bc 

Bar 

12.25 

3.8 

3.47 

1.00 

1239 

0.08 

NGC 3728 

DJ 

SAb 

Unb 

13.80 

4.2 

3.7 

0.75 

6904 

0.48 

NGC 3831 

EN 

SAB0+ 

Box* 

14.5 



0.24 

4715 

0.33 

NGC 3883 

DJ 

SAb 

Bar 

13.10 

4.1 

3.44 

0.91 

6937 

0.48 

NGC 3945 


SB0+ 

Bar 

11.38 




1220 

0.09 

NGC 4472 

EN 

E2 


9.30 




744 

0.05 

NGC 5020 

DJ 

SABbc 

Bar 

12.50 

3.6 

3.11 

0.85 

3284 

0.23 

NGC 5326 

PB 

SAa 

Unb 

12.92 

4.05 

3.97 

0.50 

2573 

0.18 

NGC 5362 

PB 

SAb 

Unb 

13.14 

3.56 

3.24 

0.37 

2314 

0.16 

NGC 5375 

DJ 

SBab 

Bar 

12.40 

3.9 

3.47 

0.81 

2418 

0.17 

NGC 5389 

PB 

SABO/a 

Box* 

13.10 

4.12 

4.10 

0.20 

1996 

0.14 

NGC 5422 

PB 

SAO 

Ell* 

12.81 

4.17 

4.09 

0.20 

1921 

0.13 

NGC 5577 

PB 

SAbc 

Ell* 

13.05 

3.84 

3.54 

0.28 

1702 

0.12 

NGC 5689 

PB 

SBO/a 

Box?* 

12.54 

4.14 

4.12 

0.25 

2295 

0.16 

NGC 5707 

PB 

SAab 

Ell* 

13.38 

4.24 

3.92 

0.25 

2354 

0.16 

NGC 5719 

PB 

SABab 

Box?* 

13.1 

4.54 

3.84 

0.36 

1676 

0.12 

NGC 5746 

PB 

SABb 

Pea* 

11.38 

4.42 

4.50 

0.16 

1676 

0.12 

NGC 5793 

EN 

SABb 

Box? 

14.30 



0.37 

3387 

0.23 

NGC 5838 

PB 

SAO- 

Box?* 

11.74 

4.21 

4.11 

0.35 

1338 

0.09 

NGC 5987 

PB 

SAb 

Ell 

13.00 

4.46 

4.14 

0.40 

3207 

0.22 

NGC 6246A 

DJ 

SABc 

Unb 

14.10 

3.9 

3.23 

0.91 

5495 

0.38 

NGC 6368 

PB 

SAb 

Ell* 

13.10 

4.84 

4.58 

0.20 

2904 

0.20 

NGC 7311 

PB 

SAab 

Ell 

13.36 

4.35 

4.07 

0.50 

4762 

0.33 

NGC 7332 

PB 

SABO 

Box* 

12.11 

3.75 

3.58 

0.26 

1584 

0.11 

NGC 7457 

PB 

SABO- 

Box? 

11.86 

3.69 

3.50 

0.52 

1115 

0.08 

NGC 7537 

PB 

SAbc 

Ell* 

13.65 

3.88 

3.62 

0.34 

2888 

0.20 
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Table 2. Spectrograph specifications during our observing runs 


Detector 

Obs. Dates (M/D/Y) 

Grating 

Disp. (A/pix) 

Approx. Resol. (A) 

Scale (Arcsec/pix) 

DIS I Blue 

1/10/00-2/11/02 

Med 

3.18 

5.7 


1.086 

DIS I Red 

1/10/00-2/11/02 

Med 

3.53 

8.6 


0.605 

DIS II Blue 

4/13/02-10/09/02 

Low 

3.05 

8.6 


0.600 

DIS II Red 

4/13/02-04/07/03 

Med 

3.13 

7.8 


0.605 

DIS III Blue 

03/06/03-02/15/04 

Low 

2.42 

7.7 


0.419 

DIS III Red 

05/29/03-02/15/04 

Med 

2.31 

6.9 


0.396 
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Table 3. Spectroscopic Observations 


Galaxy Axis PA Date (M/D/Y) Exp. Time (Sec) 

IC 267 Bar -25 


IC 302 

Bar 

8 

IC 1029 

Maj 

152 

NGC 266 

Bar 

0 

NGC 765 


15 

NGC 1642 


0 

NGC 2487 

Bar 

45 

NGC 2599 


-90 

NGC 2775 

Maj 

66 


Min 

156 

NGC 2787 

Maj 

109 

NGC 2916 

Min 

-80 

NGC 3384 

Maj 

50 

NGC 3544 

Maj 

-84 


Min 

6 

NGC 3681 

Bar 

-25 

NGC 3728 

Maj 

20 

NGC 3831 

Maj 

24 


Min 

114 

NGC 3883 

Maj 

-14 

NGC 3945 

Maj 

-22 

NGC 4472 


67 

NGC 5020 

Bar 

38 

NGC 5326 

Maj 

-44 


Min 

-134 

NGC 5362 

Maj 

-92 

NGC 5375 

Bar 

-10 

NGC 5389 

Maj 

3 


Min 

-87 

NGC 5422 

Maj 

-26 


Min 

64 

NGC 5577 

Maj 

56 

NGC 5689 

Maj 

-93 


Min 

0 

NGC 5707 

Maj 

39 

NGC 5719 

Maj 

-90 

NGC 5746 

Maj 

-9 


Min 

-99 

NGC 5793 

Maj 

-35 


Min 

55 

NGC 5838 

Maj 

42 

NGC 5987 

Maj 

-109 

NGC 6246A 


-90 

NGC 6368 

Maj 

47 

NGC 7311 

Maj 

24 



15 

NGC 7332 

Maj 

-24 

NGC 7457 

Maj 

-38 

NGC 7537 

Maj 

-100 


12/22/03 

1x2400 

10/9/02 

1x1230 

2x2400 

5/30/03 

3x2400 

9/17/02 

2x2400 

12/22/03 

3x2400 

12/1/03 

2x2400 

2/11/02 

3x2400 

2/11/02 

3x2400 

1/10/00 

3x1200 

1/10/00 

3x1200 

2/15/04 

2x2400 

12/1/03 

2x2400 

2/15/04 

1x900 

2x2400 

1/10/00 

3x1200 

4/25/00 

3x1200 

2/11/02 

3x2400 

3/6/03 

3x2400 

4/25/00 

3x1200 

5/3/00 

3x1800 

3/7/03 

3x2400 

2/15/04 

1x1200 

2x2400 

1/10/00 

1x1800 

3x1200 

3/6/03 

2x2400 

5/4/00 

4x1800 

2/11/02 

2x2400 

6/16/01 

4x2400 

4/7/03 

1x2400 

5/2/00 

1x1200 

4x1800 

5/2/00 

3x1800 

5/2/00 

2x1200 

5/3/00 

2x1800 

1/10/00 

2x1500 

3x1200 

6/17/01 

2x1500 

3x2400 

4/13/02 

2x2400 

5/29/03 

2x2400 

5/30/03 

1x1200 

3x2400 

4/17/02 

3x2400 

4/17/02 

3x2400 

5/4/00 

1x1800 

6/8/02 

2x1500 

2x2400 

5/30/03 

1x2400 

2/19/01 

1x2700 

2x1800 

2/20/01 

2x1800 

6/29/03 

1x2400 

7/1/03 

1x2100 

1x2400 

10/12/01 

1x2400 

7/3/00 

4x1800 

1x2400 

10/12/01 

1x1257 

4x2400 






